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ABSTRACT 



Cloud contraction and infall are the fundamental processes of star formation. While "blue-skewed" line profiles observed in high- 
mass star forming regions are commonly taken as evidence of infall by an ever increasing number of studies, their interpretation offers 
many pitfalls. Detecting infall via redshifted absorption in front of continuum sources is a much more direct and reliable method but 
so far mostly restricted toward absorption in the centimeter toward strong HII regions. Here we present a novel approach by probing 
absorption of rotational ammonia transitions in front of the strong dust emission of massive star-forming regions. A carefully selected 
sample of three regions with different stages of evolution is selected to study infall through the evolution of massive star-forming 
clumps. Redshifted absorption is detected toward all three sources and infall rates between 3 - 10 x 1CT 3 M Q yr~' are derived. 
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1. Introduction 

The earliest phases of massive star formation (MSF) are still 
poorly understood. We know that massive stars are born in dense 
clumps within giant molecular cloud complexes. Ultracompact 
HII regions (UCHIIRs) embedded within these clumps represent 
a key phase in the early lives of massive stars (see review 
bv lHoareir2005l) . In their environments, often hot (> 100 K), 
compact (<0.1 pc), and dense cores are found, some of 
which are believed to be in a stage prior to the formation of 
UCHIIRs dKurtz et alJ l2000t) . More recently, so-called high- 
mass protostellar objects (HMPOs) or massive young stellar 
objects (MYSOs) were reco gnized to likely re present an even 
earlier stage of MSF (e.g. iBeuther et al. 2002). Very recently, 
clumps within InfraRed Dark Clouds (IRDCs) were found to be 
promising candidat es for even earlier s tages in the formation of 
massive stars (see iMenten et al.l 12005. and references therein), 
sometimes called pre-protocluster clumps. 

One of the important results of studies of low-mass star 
formation has been the observation of infall motions (e.g. 
Bello che et al J 120021) . which give direct evidence of accretion. 
Models of high mass cluster formation also predict a large 
scale contraction as the protocluster evolves and molecular gas 
is funneled fr om the outer regions of the core into the c enter 
of the cluster dBonnell & Bate! 120021: iBonnellet all 12004). But 
toward high-mass star-for ming cores, the obser vational evidence 
of infall is still scarce (e.g. lWu & Evansll2003l) . 

Most studies rely on the observations of so-called "blue- 
skewed" line profiles: with a negative excitation gradient in a 
molecular clump, molecular lines with suitable optical depths 
(e.g. low / rotational lines of HCO + ) will show a central self- 
absorption dip whose blue peak is stronger than the red peak. 
But there are many pitfalls in the interpretation of this signature 
(see discussion in lEvansi 120031) . Especially at distances (few 
kpc) typical of MSF sites, other kinematical properties such as 
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rotation and outflow activity might easily be mistaken for infall. 
Furthermore, abundance variations in the clumps can change the 
self-absorption profile tremendously. Still, this method led to 
evidence for infalling envelopes of a quite a few of hot cores 
(IWvrowski et al.ll2006l) . 

Alternatively, absorption lines toward compact and strong 
continuum sources can be observed. In particular, centimeter 
ammonia inversion lines have been succ essfully detected in 
absorption toward bright UCHIIRs, e.g. bv lSollins et al.l (T2005) 
in G10. 6-0.4. They derive gas velocities consistent with material 
spiraling toward the c entral newly-formed stars. Similar results 
have been reported by Zhang & Ho| (1 19971) in the W51 massive 
star-forming region (see their Fig. 6), and by iBeltran et al.l 
(2006) toward G24.78+0.08. These results strongly indicate the 
presence of infall, as predicted by the accretion models. But this 
method is restricted to a fairly late state of MSF when an UCHII 
has already developed and can be used as a background source. 

At THz frequencies new opportunities arise now: Their 
high dust column densities make MSF clumps very bright 
far-infrared (FIR) continuum sources, especially in their early 
stages. Therefore, molecular lines with lower excitation than the 
(background) continuum temperature will be in absorption and 
will probe the kinematics of the clump on the line of sight toward 
the deeply embedded forming massive clusters. This signature is 
then a much more direct probe of infall and can be used through 
all embedded stages of MSF. 

Of particular interest is the ammonia molecule NH3 which 
has been extensively studied, mostly via its centimeter inversion 
lines, and prove n to be one of the best-known inters tellar 
"thermometers" (iHo & Townesl 119831: iMaret et al.l 120091) . In 
the THz range rotational transitions of ammonia can now 
be observed, in particular transitions from non-metastable to 
metastable levels, which will have low excitation temperatures 
and therefore are likely to be observed in absorption. In addition, 
ammonia is known to be one of the molecules that is the least 
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likely to fre eze out in the earliest, coldest stages of molecular 
clumps (e.g. Bergin & Langen ll997l) . 

We therefore present here a first small survey of sources 
in a range of evolutionary stages that we observed with the 
Stratospheric Observatory for Infrared Astronomy (SOFIA) in 
ammonia to measure their infall signatures. This observational 
evidence of accretion in the clumps harboring the forming 
massive (proto) clusters is essential for our understanding of the 
MSF process. 

2. Observations and data reduction 



We used the GREAlQ instrument dHevminck et all 1201 2l this 
volume) onboard SOFIA to observe the NH3 32+ - 22- line at 
1810.379971 GHz (taken from the JPL line catalog. lpickett et al.l 
1998). This line has the same lower level as the well-studied 
centimeter inversion line at 23.7 GHz of the (J, K)-(2,2) level, 
hence we know that this level is sufficiently populated to give 
rise to strong absorption signatures. 

To reach this frequency, a slight technical change was 
necessary: the modularity of GREAT allows the implementation 
of new technical opportunities between flights. A spare local 
oscillator (LO) source for the high-frequency channel allowed 
tuning to slightly lower than advertised frequencies (1810 
instead of 1820 GHz with the nominal LO). With this setup, the 
NH3 line became accessible, but only if tuned into the lower 
sideband of the receiver. GREAT operates in double sideband 
mode and parts of the image sideband then cover a region of very 
low transmission, which led to uncertainties in the determination 
of the (weak) continuum level. The lower frequency channel was 
tuned to the CO (13-12) line at 1496.922909 GHz and was in 
turn used, at much better transmission, to determine the THz 
continuum flux of the sources. 

The beam sizes at the observing frequencies are 20 and 16 
arcsec, respectively. The three sources given in Table Q~| were 
observed on 2011 July 20 at an altitude ranging from 39000 to 
43000 feet. 

The system temperatures during the observations and 
integration times are given in Table Q] As backend s, MPIfR fast 
Fourier transform spectrometers dKlein et al.ll2012l) were used to 
cover a bandwidth of 1.5 GHz for each channel with a velocity 
resolution of 0.03 km/s, which was later lowered to 1 km/s to 
increase the signal-to-noise ratio in the spectra. The wobbling 
secondary was chopped with a throw of 120" about the cross 
elevation axis in a symmetric mode. The pointing was monitored 
with the fiber optic gyroscopes. Spectra were calibrated to a 
T* h scale and then the conversion to Jmb was made assuming 
a forward efficiency of 0.95 and beam efficiencies of 0.54 and 
0.5 1 for the lower and higher frequency receivers, respectively. 

In addition, complementary ground-based spectra of the 
C I7 (3-2) and HNC (4-3) were obtained with the APEX 
telescope. 

Final processing and analysis of the data was then performed 
using the CLASS program within the GILD AS software Q. 

3. Results 

The averaged and baseline-subtracted spectra are shown in 
Fig. Q] Ammonia absorption lines were detected in all three 
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1 GREAT is a development by the MPI fiir Radioastronomie and 
the KOSMA/Universitat zu Koln, in cooperation with the MPI fiir 
Sonnensystemforschung and the DLR Institut fiir Planetenforschung. 

2 http://www.iram.fr/IRAMFR/GILDAS 



Fig. 1. NH3 32+ - 22- spectra of the observed sources. Results 
of Gaussian fits to the line are overlaid in green. The systemic 
velocities of the sources, determined using C n O (3-2), are 
shown with dotted lines. 



sources. The line parameters from Gaussian fits are given in 
Table [2] Compared with C lv O (3-2) line velocities, which 
as optically thin emission lines probe the systemic velocities 
of the sources with a comparable resolution of 18", all 
ammonia absorption lines are redshifted. This is consistent 
with blue-skewed emission profiles observed in HNC (4-3) 
toward these clumps with the APEX telescope. The fact that 
we observe absorption lines is not surprising: compared to 
the microwave inversion lines, the FIR rotational lines have 
Einstein A values about five magnitudes higher, hence they 
will decay rapidly from the non-metastable (J>K) into the 
metastable (J=K) levels, which in turn only decay slowly via 
collisions dHo & Towneslfl983l) . Thus, the metastable levels are 
"overpopulated" compared with the non-metastable ones. 

Because we were unable to measure the 1.8 THz continuum 
reliably owing to the double sideband reception and the low 
atmospheric transmission in part of the image band, the 1 .5 THz 
continuum was used to estimate the continuum level at the NH3 
frequency. Table[3]gives the 1 .5 THz continuum level, which was 
divided by 2 to correct for the single sideband calibration that 
by default is applied to the lines. The modeling of the clumps 
described in the next section suggests that the real 1.8 THz 
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Source 


Stage 


R.A.(J2000) Decl.(J2000) 
( h m s ) ( " " ) 


Lbol 


d 
(kpc) 


S(870 y um)° 
(mJy/bm) 


T 

(K) 


hnt 

(min.) 


W43-MM1 

G31.41+0.31 

G34.26+0.15 


MIR-quiet core 
hot core 
UCHII 


18 47 47.0 -0154 28 
18 47 34.3 -0112 46 
18 53 18.6 +0114 58 


2.3 x 10 46 

2.6 x 10 5c 

4.7 x 10 5c 


5.5 fc 
1.9 C 
3.7 C 


21.2 
21.2 
44.7 


3700 
3500 
3500 


13 
6.8 
2.5 


















ISchuller et al.l d2009h. iMotte et al 


.1 i2003h. iRolffs et alJ d201 lh. 





Table 2. Line parameters from Gaussian fits to the NH3 lines. 
Nominal fit errors are given in parenthesis. In addition, the 
velocity of C 17 (3-2) lines observed with the APEX telescope 
are given. 



Source 



(K) 



Av 
(kms- 1 ) 



JW 

V LSR 

(km s~ 



„c"o 

V LSR 

(kms^ 1 ) 



W43-MM1 

G31.41+0.31 

G34.26+0.15 



-0.96 (0.22) 
-1.18(0.29) 
-3.38 (0.56) 



5.3 (0.8) 
3.7 (0.8) 
5.5 (0.6) 



99.7 (0.4) 
99.4 (0.4) 
61.2(0.3) 



97.65 (0.06) 
97.02 (0.04) 
58.12(0.03) 



Source 



-* cont 

(K) 



N(2, 2) 
(10 14 cnr 2 ) 



W43-MM1 

G31.41+0.31 

G34.26+0.15 



1.5 
1.3 

5.5 



1.0 

2.4 
1.0 



1.2 
2.0 
1.2 



continuum might be up to 20% higher. The optical depths of 
the absorption can then be estimated from the line-to-continuum 
ratio as r = -log(l + Tmb/Tcom) using the temperatures from 
Table [2] and [3] Optical depths on the order of unity are found 
meaning that the filling factor of the absorbing material in front 
of the dust continuum must be very high. 

The total colum n density of th e absorbing NH3 can then be 
computed from (e.g.|Comito et al. 2003b, 



N 2 - = 

2 



8ttv 3 gi_ 



— tAv. 



(1) 



This equation is valid for excitation temperatures much 
lower than the continuum temperature and also much lower 
than the 'transition temperature' T u \ (hv/k — 86 K for this 
NH3 transition), which likely is the case given the rapid 
radiative decay of the non-metastable upper level. The transition 
probability of A„i - 3.56 x 10~ 2 s leads to a critical density for 
the transition of 10 9 ctrT 3 . gi and g„ are the statistical weights 
of the lower and upper levels, respectively. Both inversion levels 
have almost the same population for densities above 10 4 cm -3 , 
hence the total column density in the (2,2) levels is twice as high. 
This leads to column densities of 1 - 2 x 10 14 ctrr 2 as given in 
Table [3] 

In Figure|2]we sh ow a comparison of the (2,2) in version line 
seen with the VLA dKeto et all 1 19871 iGomez et alJl2000h with 




3 2+ -2 2 _ ® 1810 GHz 

2 2 .,-2 2 _ « 23 GHz (Gomez+2000; VLA) 



Table 3. Continuum temperatures at 1 .5 THz, the NH3 optical 
depths and the column densities in the (2,2) levels. 
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Fig. 2. G34.26+0.15 SOFIA NH 3 spectrum compared with 
the VLA NH3 (2, 2) spectrum (shown in red) taken from 
IGomez et al. (2000), which was integrated over the region that 
shows absorption. A two-component hyperfine fit to the (2,2) 
spectrum is shown in green. 



the SOFIA results. This inversion line shares the same lower 
level with the line observed by SOFIA. The velocity and width 
of the two absorption lines agree very well. Interestingly, the 
VLA spectrum shows the (2,2) line's outer hyperfine structure 
components in emission. This can be understood when the 
absorbing continuum is only filling the beam partly and the 
absorption spectrum is blended with an additional high optical 
depth emission component from the hot core that is situated 
offset from the continuum. 



4. Analysis 

To model the observed absorptio n line s , we u sed the physical 
structures determined by IRolffs et all (1201 lh : they fitted the 
continuum structure of MSF clumps as observed in the 
ATLASGAL survey (ISchuller et al.l 120091) with density power 
laws n ~ r~ L5 for our observed sources) and determined 
temperature structures dictated by the inner heating sources. 
The velocity structure was assumed to be a fixed fraction 
of the free-fall velocities and adjusted using RATRAN 
(Hogerheijde & van der Tak 2000) to fit high-density molecular 
probes observed with the APEX telescope. We compared the 
predicted 1 .5 THz continuum of the mo dels with the observed 
ones. G3 1.41 was fitted well with dust of Ossenkopf & Henning 
(1994) using their models with coagulation at n - 10 5 cm -3 
but no ice mantles. For G34.26 we adjusted the grain properties 
to dust with thin mantles that have a higher spectral index 
to reproduce the high obser ved continuum tem perature. W43- 
MM1 was not modeled by Rolf fs et all (1201 lh . For this source 



Wyrowski et al.: Terahertz ammonia absorption as a probe of infall in high-mass star forming clumps 



we scaled the G31.41 model to match the observed 1.5 THz 
continuum flux. 

In the next step we adjusted the ammonia abundance 
(n(H 2 )/n(NH 3 )) and the velocity to fit the SOFIA NH 3 
absorption lines with the RATRAN code. While in principle 
the clumps could have an increased hot core NH3 abundance 
in the inner T > 100 K region, this hot core component to the 
spectrum is negligible because of its small filling factor. To fit 
the velocity, we adjusted the turbulence widths of the gas and 
the free-fall fraction / i n v — — / X V2 GMj n /r, where M- m 
is the mass inside r ("cf. iRolffs et alJl201lh . An example fit is 
shown in Fig.[3]forG34.26+0.15 together with a fit of the model 
to the NH3 (2,2) inversion line observed by Churchwell et al.l 
(1990J). In the model, the velocity was adjusted to 30% of the 
free-fall velocity and an ammonia abundance of 0.8 x 10~ 8 . The 
Doppler ^-parameter was adjusted to b - 1.5 km s _1 . In the 
SOFIA spectrum, there is even evidence for additional excess 
redshifted absorption, not fit by the model, albeit only with low 
signal-to-noise ratio. The optically thick main cm-inversion line 
is somewhat overpredicted, whereas the outer hyperfine satellites 
are reproduced quite well. High infall velocities in the innermost 
part of the model do not contribute to the observed line profile, 
since the filling factor of the hot core region is too low. W43- 
MM1 and G31.41 could be fitted with velocity fields with 20% 
free-fall contribution. 




50 60 70 

Velocity (km/s) 

Fig. 3. G34.26+0 15 N H 3 SOFIA and Effelsberg (taken from 
IChurchwell et al.1 Il990h spectra in red compared with the 
RATRAN model discussed in the text (black profiles). The 
systemic velocity from C I7 is indicated by a dotted line. 



5. Discussion and conclusions 

From the observed and modeled velocity fields the infall rate 
of the clumps can be determined. From the RATRAN models, 
M was determined from the model density and velocity at the 
radius corr esponding to the SOFIA beam size. Alternatively, 
Eq. 4 from iBeltran et al. (2006) can be used with the observed 
column density and infall velocity. Both approaches agree and 
lead to infall rates between 3 - 10 x 10~ 3 Mp yr' 1 . These rates 
are high enough to quench any HII region dWalmslevI [1995) 



and to overcome the radiation pressure from the forming (pro to) 

cluster, although we only probed the infall at the outer envelope. 

In the inner part accretion might already be halted. From the 

absorption, we also only probed the infall on the line-of-sight, 

hence it might not be spherical. 

In summary, we present the first velocity-resolved ammonia 

spectra from a rotational transition at THz frequency toward 

MSF clumps. Absorption of a comparable NH3 3 2 - - 2 2+ line 

at 1764 GHz was also recently obser ved with Herschel in the 

low-mass source IRAS 16293-2422 dHilv-Blant et alJl2010h . In 

G34.26 the redshifted absorption agrees with cm absorption 

observed with the VLA toward the UCHII component B 
■ 11 —1 

(Heaton et al. 1989). This component is within the SOFIA beam 
but the strong absorption we detected implies that the absorption 
fills most of the dust continuum, which peaks toward the hot 
molecular core. G31.41 and W43-MM1 do not have compact 
cm continuum that might be used for absorption studies at 
cm wavelengths, hence our observations are first detections 
of redshifted absorption in these sources and demonstrate that 
ammonia absorption in front of strong dust continuum can be 
successfully used to probe infall in a variety of evolutionary 
stages of MSF clumps. 

Acknowledgements . Based [in part] on observations made with the NASA/DLR 
Stratospheric Observatory for Infrared Astronomy. SOFIA Science Mission 
Operations are conducted jointly by the Universities Space Research 
Association, Inc., under NASA contract NAS2-97001, and the Deutsches SOFIA 
Institut under DLR contract 50 OK 0901. 



References 

Belloche, A., Andre. P., Despois, D., & Blinder, S. 2002, A&A, 393, 927 

Beltran, M. T., Cesaroni, R., Codella, C, et al. 2006, Nature, 443, 427 

Bergin, E. A. & Langer, W. D. 1997, ApJ, 486, 316 

Beuther, H., Schilke, P., Menten, K. M., et al. 2002, ApJ, 566, 945 

Bonnell, I. A. & Bate, M. R. 2002, MNRAS, 336, 659 

Bonnell, I. A., Vine, S. G., & Bate, M. R. 2004, MNRAS, 349, 735 

Churchwell, E., Walmsley, C. M., & Cesaroni, R. 1990, A&AS, 83, 119 

Comito, C, Schilke, P., Gerin, M., et al. 2003, A&A, 402, 635 

Evans, II, N. 2003, in SFChem 2002: Chemistry as a Diagnostic of Star 

Formation, ed. C. L. Curry & M. Fich, 157 — h 
Gomez, Y., Rodriguez-Rico, C. A., Rodriguez, L. F., & Garay, G. 2000, Rev. 

Mexicana Astron. Astrofis., 36, 161 
Heaton, B. D., Little, L. T., & Bishop, I. S. 1989, A&A, 213, 148 
Heyminck, S., Graf, U., Glisten, R., et al. 2012, this volume 
Hily-Blant, P., Maret, S., Bacmann, A., et al. 2010, A&A, 521, L52+ 
Ho, P. T. P. & Townes, C. H. 1983, ARA&A, 21, 239 
Hoare, M. G. 2005, Ap&SS, 295, 203 

Hogerheijde, M. R. & van der Tak, F. F. S. 2000, A&A, 362, 697 
Keto, E. R., Ho, P. T. P., & Reid, M. J. 1987, ApJ, 323, LI 17 
Klein, B., Hochgiirtel, S., Kramer, I., Bell, A., & Meyer, K. Giisten, R. 2012, 

this volume 
Kurtz, S., Cesaroni, R., Churchwell, E., Hofner, P., & Walmsley, C. M. 2000, 

Protostars and Planets IV, 299 
Maret, S., Faure, A., Scifoni, E., & Wiesenfeld, L. 2009, MNRAS, 399, 425 
Menten, K. M., Pillai, T., & Wyrowski, F. 2005, in IAU Symposium, Vol. 227, 

Massive Star Birth: A Crossroads of Astrophysics, ed. R. Cesaroni, M. Felli, 

E. Churchwell, & M. Walmsley, 23-34 
Motte, F, Schilke, P., & Lis, D. C. 2003, ApJ, 582, 277 
Ossenkopf, V & Henning, T. 1994, A&A, 291, 943 
Pickett, H., 

Poynter, R. L., Cohen, E. A., et al. 1998, J. Quant. Spec. Radiat. Transf., 

60, 883 
Rolffs, R., Schilke, P., Wyrowski, E, et al. 201 1, A&A, 527, A68+ 
Schuller, F, Menten, K. M., Contreras, Y, et al. 2009, A&A, 504, 415 
Sollins, P. K., Zhang, Q., Keto, E., & Ho, P. T. P. 2005, ApJ, 624, L49 
Walmsley, M. 1995, in Revista Mexicana de Astronomia y Astrofisica, vol. 27, 

Vol. 1 , Revista Mexicana de Astronomia y Astrofisica Conference Series, ed. 

S. Lizano & J. M. Torrefies, 137 
Wu, J. & Evans, II, N. J. 2003, ApJ, 592, L79 

Wyrowski, E, Heyminck, S., Giisten, R., & Menten, K. M. 2006, A&A, 454, L95 
Zhang, Q. & Ho, P. T. P. 1997, ApJ, 488, 241 



